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Abstract 

The first observation is made of hadronic string breaking due to dynamical 
fermions in zero temperature lattice QCD. The simulations are done for SU(2) 
color in three dimensions, with two flavors of staggered fermions. The results 
have clear implications for the large scale simulations that are being done to 
search (so far, without success) for string breaking in four-dimensional QCD. 
In particular, string breaking is readily observed using only Wilson loops to 
excite a static quark-antiquark pair. Improved actions on coarse lattices are 
used, providing an extremely efficient means to access the quark separations 
and propagation times at which string breaking occurs. 
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The string picture of quark confinement predates Quantum Chromodynamics, and con- 
tinues to play a central role in theoretical efforts to characterize the physics of QCD [|TJ. 
Simulations of quenched lattice QCD have in fact demonstrated that color-electric field 
lines connecting a static quark and antiquark are squeezed into a narrow tube [fj] , in accord 
with many models of the linearly rising static potential. 

In QCD with dynamical fermions the flux-tube joining static quarks is expected to be 
unstable against fission at large R, where there is enough energy in the fields to materialize 
light quarks from the vacuum, which bind to the heavy quarks to form a pair of color- 
neutral bound states. The simplest indicator of string breaking is provided by the static 
quark- antiquark potential V(R), which should approach a constant at large R 

V(R - oo) = 2M Q5 , (1) 

where Mqq is the mass of a bound state of a heavy quark and a light antiquark. 

There is as yet no convincing evidence of string breaking in lattice simulations at zero 
temperature, despite extensive simulations by several large scale collaborations 0^-^]. A 
recent conjecture f7||l, which has received widespread attention, holds that the failure to 
detect the broken string state is due to a very poor overlap of the Wilson loop operator 
with the true ground state of the system at large R (Wilson loops have been used almost 
exclusively to construct trial states in these simulations). The possibility has been raised that 
the broken string state may even be unobservable with Wilson loops for all practical purposes 
0. This picture has recently stimulated interest in simulations using other operators to 
excite the system P,|lO|. 

Evidence for a somewhat different point of view is presented in this paper. The task of 
observing string breaking in full QCD is computationally very challenging, given the high 
cost of simulating dynamical fermions. These simulations have so far been done on lattices 
with a relatively fine grid. Consequently, only a fairly limited range of quark separations 
R has been explored and, more importantly, the propagation of the trial states created by 
Wilson loops has only been studied for rather short Euclidean times T (e.g., T = 0.4-0.8 fm 
in Ref. [[J). The restriction to such short propagation times precludes a definitive assessment 
of whether string breaking occurs in these states. This limitation may be substantially 
alleviated by working on coarser lattices with improved actions. 

This proposal is assessed here by doing simulations in three-dimensional lattice QCD 
(QCD3) with two flavors of dynamical staggered fermions. QCD3 provides an excellent 
laboratory for studying the physics of realistic QCD (for a review see Ref. [PH ). Quenched 
QCD3 has been shown to exhibit linear confinement, flux-tube formation, a deconfining 
phase transition, and a rich glueball spectrum, and spontaneous breakdown of a "chiral" 



symmetry has been observed in both the quenched and unquenched theories . 

The first observation of string breaking by dynamical fermions is reported here, using 
only Wilson loop operators to excite the static quark-antiquark system |R|. Although 



these simulations are done in QCD3, the results have clear implications for work on the 
four- dimensional theory. In particular the simulations are done on coarse lattices, using an 
improved gluon action that is accurate through 0(a 2 ), and a staggered fermion action that 
is accurate through 0(a). This allows the system to be much more efficiently probed at 
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the physical separations and propagation times relevant to string breaking. Moreover, the 
overlap of an operator with the lowest-lying state generally increases with the coarseness of 
the lattice, due to the suppression of higher momentum modes in the trial state. In contrast, 
recent studies of QCD 3 with dynamical scalar matter fields ||, and adjoint matter fields 



I4| , using unimproved actions on lattices with about half the spacing considered here, did 
not resolve string breaking in Wilson loops, presumably because insufficient propagation 
times were attained. 

Using the quenched string tension to set the "physical" length scale for QCD 3 , these 
results suggest that string breaking occurs in Wilson loops at modest propagation times, 
somewhat in excess of 1 fm (on modestly coarse lattices), which is not much beyond what 
has been reached in four dimensions. The physical quark separation R « 1.5 fm at which 



string breaking is found in QCD3 is consistent with expectations in four dimensions ||T5 
Lattice spacings of about 0.2-0.3 fm are used here; these lattices are about twice as coarse as 
those that have been used in four- dimensional studies of string breaking. Comparable lattice 
spacings could be reliably accessed in realistic QCD using improved actions |16|], resulting 
in an enormous improvement in computational efficiency 

To begin with, consider the pure gauge theory in three dimensions. To reduce the 
computational cost SU(2) color is considered. Simulations were done using a tree-level 
0(a 2 )-accurate gluon action, allowing for different spatial and temporal lattices spacings, a s 
and a t respectively fl~7| 



^imp (3 ^ t £/j,u 3-^^ 12 (^l tl/ -^W*) ' (^) 

P^ v is one-half the trace of the lxl plaquette and R^ u is one-half the trace of the 2x1 
rectangle in the fi x v plane. The bare QCD coupling g\ (of dimension 1) enters through 
the dimensionless parameter (3 = A/g^ag, and the bare lattice anisotropy is input through 
&i = &3 = a g /at and = a t /a s = 1, 2). 

QCD3 is a super- renormalizable theory; in fact the bare coupling g\ and bare quark 
masses mo remain finite as the ultraviolet regulator is removed. Since g^ and mo are cutoff- 
independent in the continuum limit, one obtains extremely simple scaling laws for physical 
quantities in lattice QCD 3 . In particular, masses in lattice units (including the input bare 
quark masses) should satisfy 

(3 am = constant (3) 

in the continuum limit. The improved action Si mp differs from the continuum theory at the 
tree-level by terms of 0(a 4 ); one- loop corrections induce errors of O(go<2 3 ). 

Results for the quenched string tension obtained with S imp are shown in Fig. [1]; for 
comparison, results obtained by Teper [|ITJ for the unimproved Wilson gluon action are also 



shown (the dashed lines show the continuum limit estimated in ||11|| ). A comparison of the 
potential computed from the two actions on lattices of comparable spacing is shown in Fig. 

B 

Simulations of S- imp were done on lattices with /3 = 2, 2.5, and 3, all with anisotropy 
a t/°s = 1/4- The lattice volumes range from 16 3 at (3 = 2 to 24 3 at f3 = 3. Hybrid 
molecular dynamics were used to generate the configurations (the $ algorithm [18j was used 
to generate the unquenched configurations analyzed below); 50 molecular dynamics steps 
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were taken with step size At = 0.02. A standard fuzzing procedure was employed on the 
link variables used in Wilson loop measurements. Integrated autocorrelation times satisfied 
r int ~ 0.5, with typically 5-10 trajectories skipped between measurements. 

A superficially surprising aspect of these results is the 0(a) scaling violation evident in 
the Wilson action data for \fa j g 2 . In fact this demonstrates the need to renormalize the bare 
coupling at finite a [plj . In four dimensions it is known that perturbative expansions in the 
bare lattice coupling are spoiled by large renormalizations. A renormalized coupling defined 
by a physical quantity absorbs these large corrections, and greatly improves perturbation 
theory for many quantities [20|. 



In lattice QCD 3 one expects to find an 0(<7qo) renormalization of the bare coupling at 
one-loop order. This produces the linear scaling violation in the unimproved results for 
y/a jg\. On the other hand this effect can be removed by using a physical quantity to define 
a renormalized scale. For example, Teper has done simulations of three-dimensional glueball 
masses mc and showed that the ratios m^/ ' \fo exhibit 0(a 2 ) scaling violations. 



The distinctive signature of the renormalization of the bare coupling in QCD3 exposes an 
interesting feature of the action. One might expect to see a reduction in the renormalization 
of (?q when S imp is used, compared to simulations with the Wilson action. Remarkably one 
finds that the renormalization is in fact almost completely eliminated; as seen in Fig. [1| the 
error in \fo jg 1 at ft = 2 is reduced from about 45% with the Wilson action, to less than 
about 5% with Si mp . This is also apparent from the data in Fig. |[ where the potential 
V and separation R from the lattices at ft = 2 and 3 are compared on a common scale, 
set by (/q. This is a genuinely surprising result, since all of the higher dimension operators 
present in Si mp should contribute to a leading O(g^a) renormalization of the bare coupling. 
Apparently the operator series in this effective action converges very rapidly, even at scales 
near the lattice cutoff. 

One can therefore conveniently measure physical quantities directly in terms of the bare 
coupling when the improved action Eq. (^) is used. It is important to note however that 
the use of S imp would prove advantageous even if g\ required renormalization; one would 
simply use a measured quantity (such as the quenched string tension) to set the scale for 
other observables, and the leading 0(a 2 ) errors in the Wilson action would thus be removed. 

Simulations with dynamical Kogut-Susskind fermions were done to look for effects of 
string breaking. The staggered fermion action iSks i n three dimensions [^TJ is identical in 
form to the four- dimensional theory, and describes two flavors of four-component spinors 



Sks = V»( x )x(x) [U^(x)x(x + fx) 
-UUx - fx)x(x - /})] + 2am J2 x(x)x(x), (4) 



where r)^(x) = (— +^-1 is the usual staggered fermion phase. Unquenched simulations 
were done with improved glue on isotropic lattices (a t = a s ) at ft = 2 and 3, on 12 2 x 8 
and 16 2 x 10 volumes respectively f22 |. Isotropic lattices for the unquenched simulations 



were used in order to probe the longest propagation times T possible, for the least computa- 
tional cost (although measurements over a finer range of T on anisotropic lattices would be 
advatangeous for fitting and possible extraction of excited state energies). The input bare 
quark mass in lattice units was scaled according to m / g$ = 0.075 [cf. Eq. (B)], with m^/rrip 
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found to be m 0.75. Approximately 6,000 measurements were done at (3 = 2, and 2,000 at 
P = 3. 

The leading discretization errors in the action Eq. (H) are of 0(a 2 ) P3"| ; at finite lattice 
spacing the bare quark mass must also absorb an O^a) renormalization. However these 
effects appear to induce little scaling violation in the static potential at the small quark 
mass used here, perhaps because the energy in the color fields dominates the string breaking 
process (significant scaling violations in the unquenched potential were seen when the Wilson 
gluon action was used). 

Results for the unquenched potential are shown in Fig. [5]. The onset of string breaking 
is assessed by comparing the potential energy with the mass of the two Qq bound states into 
which the system should hadronize. The staggered heavy-light meson propagator |24| was 
computed in the unquenched configurations at (3 — 3, for equal valence and sea quark masses; 
the result is shown as the dashed lines in Fig. |3|. String breaking is clearly demonstrated, 
with the unquenched potential approaching the expected asymptotic value [cf. Eq. (]]])]. 
Excellent scaling behavior is observed; note that the lattice spacing increases by 50% from 
/3 — 3 to /3 — 2. A continuum extrapolation of the string breaking distance Rb can be 
estimated from these results, Rb/(f3a) « 2.5. 

The onset of string breaking is made particularly evident by a comparison of the quenched 
and unquenched potentials at the same j3 (note that no adjustment of the energy zero has 
been made in any of these results). The meaningfulness of this comparison is supported 
by the fact that the bare coupling g$ undergoes little renormalization; moreover quenching 
effects at a given f3 were found to change the p-meson mass by less than 10%. 

Only (fuzzy) Wilson loop operators were used in these calculations. As discussed above, 
little evidence of string breaking in Wilson loops has been found in previous simulations 
with dynamical fermions. A distinguishing feature of the simulations done here, compared 
with earlier work, is the relative coarseness of the lattices that were used (of course the 
physics of the string breaking process could be substantially different in three dimensions; 
on the other hand, simulations of QCD3 on much finer lattices |]9,|l4| did not resolve string 
breaking in Wilson loops). 

Using coarse lattices made it much easier to measure Wilson loops W(R,T) at propaga- 
tion times T large enough to accurately identify the ground state energy at large R. Plots 
of the time-dependent effective potential V(R,T) = —\n[W(R,T)/W(R,T — 1)] are given 
in Fig. |], and show that one could easily be misled as to the shape of the potential if the 
calculation is not done at large enough T, particularly at large R. 

In order to have some intuition for the length scales probed in these simulations, it is 
instructive to use the physical value of the string tension in four dimensions, a/ct = 0.44 GeV, 
to in effect set the value of the dimensionful coupling constant in QCD3. Using the quenched 
continuum extrapolation f3a^/a = 1.33(1) |Tl| one identifies, for example, a(/3 = 3) ~ 0.2 fm. 



Comparable, if slightly smaller, estimates of the lattice spacing are obtained if simulation 
results for the quenched or unquenched p meson mass in QCD 3 are identified with the 
physical mass. 

This scale setting implies a string breaking distance Rb ~ 1.5 fm, which is numerically 



very similar to estimates of Rb in four dimensions The results of Fig. |] suggest 

that propagation times somewhat in excess of 1 fm are needed to adequately resolve string 
breaking in Wilson loops on modestly coarse lattices; notice than an appreciable rise in the 
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potential remains even at T rs 0.6 fm, which is comparable to the longest propagation times 
resolved in the four- dimensional simulations of Ref . || . 

Although longer propagation times T are necessary to fully isolate the ground state en- 
ergy at the largest separations R in Fig. £|, these data nevertheless clearly support the onset 
of string breaking. This follows from a comparison of the unquenched and quenched poten- 
tials. Of particular importance is the fact that any slope that may exist in the unquenched 
potential data at large R is clearly much smaller than the slope in the quenched data. Such 
a small slope, if truly present, would imply a very large change in other physical quantities 
due to unquenching. In fact, m p was found to change by less than 10% with unquenching. 
Hence the significant flattening of the unquenched potential data in Figs. |3| and |], compared 
to the quenched data, is very strongly indicative of the onset of string breaking. 

To summarize, string breaking due to dynamical fermions was observed in three- 
dimensional lattice QCD. The results have clear implications for current work on the four 
dimensional theory. In particular, string breaking was readily observed using only Wilson 
loop operators. The use of improved actions on coarse lattices provided a crucial advantage 
in accessing the quark separations and propagation times at which string breaking occurs. 
The use of comparably coarse lattices in four dimensions would result in an enormous im- 
provement in computational efficiency, and could enable string breaking to be demonstrated 
in realistic QCD. 

I am indebted to Peter Lepage for several insightful discussions. I also thank R. 
Woloshyn, R. Fiebig, M. Alford, and M. Teper for fruitful conversations. This work was 
supported in part by the Natural Sciences and Engineering Research Council of Canada. 
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FIGURES 



FIG. 1. Lattice spacing dependence of the quenched string tension from the improved gluon 
action (□); results for the Wilson gluon action (•) are taken from Ref. |11]. 



FIG. 2. Static potential from the Wilson gluon action at (3 = 2 (•), and from the improved 
gluon action at (3 = 2 (□) and (3 = 3 (x). 



FIG. 3. Static potential from the improved gluon action with and without dynamical staggered 
fermions, at j3 = 2 (□) and j3 = 3 (x). The unquenched simulations were done on isotropic 
lattices (a t = a s ). The dashed lines are one sigma limits for 2(3a s MQq, computed on dynamical 
configurations at (3 = 3. 



FIG. 4. Time-dependent effective potential V(R,T) versus R, in the dynamical configurations 
at (3 = 3, for T/a s = 1 (□), 3 (o), and 5 (A). Also shown is the quenched potential at (3 = 3 (•), 
and the broken string energy (dashed lines). 
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